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We study the photoluminescence spectra of n-doped InP bulk wafers, both in the reflection and 
the transmission geometries relative to the excitation beam. From the observed spectra we estimate 
the spatial distribution of minority carriers allowing for the spectral filtering due to re-absorption 
of luminescence in the wafer. This distribution unambiguously demonstrates a non-exponential 
drop-off with distance from the excitation region. Such a behavior evidences an anomalous photon- 
assisted transport of minority carriers enhanced owing to the high quantum efficiency of emission. It 
is shown that the transport conforms very well to the so-called Levy-flights process corresponding to 
a peculiar random walk that does not reduce to diffusion. The index 7 of the Levy flights distribution 
is found to be in the range 7 = 0.64 to 0.79, depending on the doping. Thus, we propose the high- 
efficiency direct-gap semiconductors as a remarkable laboratory system for studying the anomalous 
transport. 

PACS numbers: 78.30.Fs,78.55.Cr,78.60.Lc,74.62.En 

Keywords: photoluminescence, photon recycling, anomalous diffusion 



I. INTRODUCTION 

The shape of the luminescence spectra at room and 
cryogenic temperatures is known to be strongly affected 
by the doping: it modifies the absorption spectrum and 
therefore both the emission spectra and the radiation es- 
cape. More detailed physical picture of luminescence in 
heavily doped samples involves a number of different fac- 
tors, such as change in the density of states, the bandgap 
narrowing, filling of the impurity and the conduction 
band states, carrier scattering and the resultant momen- 
tum non-conservation in optical transitions, as well as 
the effects of diffusion of the minority carriers followed 
by surface recombination.— The relative impact of these 
factors is difficult to discriminate in heavily doped sam- 
ples (Nd > 10 18 cm~ 3 ) and equally good fitting of the 
spectra can be obtained with different set of parameters 
(see Refs. [^-13 ) leading to different interpretations of the 
experimental data. 

At moderate doping levels (Nd < 10 18 cm -3 ) the 
electron system is non-degenerate and the luminescence 
spectra are not modified by the conduction band filling. 
The spectra remain similar for a whole range of moderate 
doping concentrations. The new phenomenon that comes 
into play in high-quality moderately-doped samples is 
the "photon recycling" that is the multiple radiative 
emission-reabsorption of luminescence photons^— The 
radiative recombination processes in such samples are 
several orders of magnitude faster than the non-radiative 
recombination. Besides, the residual free-electron ab- 
sorption of interband luminescent radiation becomes neg- 
ligible. With negligible losses of minority carriers or pho- 
tons, the kinetics of minority carrier luminescence is dom- 
inated by the "photon recycling" . Two main features of 
the recycling have been noted£ (i) the enhanced appar- 
ent lifetime of the minority carriers and (ii) enlarged mi- 
nority carrier spread in the sample beyond the excitation 
area. Traditionally, the enhanced spread is interpreted in 



terms of a modified photon-assisted diffusivity and stud- 
ied by numerical calculations? 8 - which do not reveal the 
underlying physics. However, as will be shown in the 
present work, the spatial distribution of minority carri- 
ers is not exponential and the diffusion approximation 
fails completely. 

In the photon-assisted diffusion model, the distribution 
of minority carriers is formed by two additive transport 
processes: (i) the random flights of holes (at sub-micron 
distances) interrupted by scattering, as in ordinary dif- 
fusion, and (ii) the photon-assisted transport of holes at 
much larger distances. The latter process corresponds 
to radiative recombination of a hole with the resultant 
interband photon, when reabsorbed, producing another 
hole at a different spot. If the rcabsorption length is 
strictly limited, e.g. by the residual non-interband ab- 
sorption, then the photon-assisted process can be viewed 
as a random walk with enlarged steps. On a large scale, 
this would be equivalent to a diffusive particle spread 
with much enhanced diffusivity. In this case, the hole 
distribution should decay exponentially away from the 
excitation area. The conventional hole diffusion mecha- 
nism implies a small diffusion length Ld < 10 /im, but 
with photon recycling one can envision at least a tenfold 
enhanced Ld- 

However, in the experimental situation studied in this 
work, the effective Ld — > 00 and the diffusion approxi- 
mation fails as a matter of principle. The free path of 
photons depends on their wavelength and even though 
the "typical" interband photons are re-absorbed at short 
distances, for those emitted in the red wing of the lu- 
minescence spectrum the free path grows exponentially 
and much exceeds the sample size. The overall distribu- 
tion of the reabsorption probability and the distribution 
of the secondary holes created by reabsorption both ac- 
quire long non-exponential tails. The resultant random 
process of the temporal spread of the minority carrier 
concentration substantially differs from the normal diffu- 
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sion and should be considered in terms of Levy flights, 9 
The general Levy-flights (LF) transport corresponds to 
a random walk with the distribution of jump lengths 
V{x) characterized by a divergent second moment, e.g. 
V{x) oc |cc|~( 1+7 ) with the exponent < 7 < 2. The con- 
ventional diffusivity is not even defined for such a random 
walk. 

The LF transport problem has been extensively stud- 
ied mathematically. Description of the anomalous trans- 
port in terms of fractional dynamical equations or, for 
random walks in an external field, fractional Fokker- 
Planck equations is amply discussed in the reviews^^ 
LF phenomena are well-known in astrophysics, as they 
occur in the problem of transport of resonance radia- 
tion in celestial bodies^ they are also known in plasma 
physics as the imprisonment of resonance radiation in 
gaseous discharge ! 12 ' 13 Interestingly, these phenomena 
are more common than one would think: thus, Levy 
flights were recently invoked to explain movement strate- 
gies in mussels as revealed in the patterning of mussel 
beds^ as well as ocean predators search strategies in 
regions where prey is sparse 

Here we report the results of experimental and the- 
oretical studies of the luminescence spectra of different 
n-type InP wafers. We demonstrate that these spectra 
are well suited to study the anomalous transport of the 
minority carriers. The spatial distribution of minority 
carriers is revealed by the analysis of the ratio of spectra 
observed in the transmission and reflection geometries. 
The observed spectra are distorted due to the energy- 
dependent re-absorption of outgoing radiation; we quan- 
tify this effect by filtering functions. Because of the ex- 
ponential dependence of the absorption coefficient near 
the bandgap, the input of remote layers to the observed 
luminescence is strongly filtered and shifted to the red 
wing of the spectrum. The energy-dependent filtering 
makes the experimentally observed spectra sensitive to 
the spatial distribution of minority carriers. 

The direct manifestation of the filtering effect is the 
red shift of the observed position of the maximum in 
the emission spectrum. This shift was clearly seen in 
the evolution of the luminescence line observed from the 
side edge of the sample when the distance from the exci- 
tation spot to the sample edge is varied up to several 
centimeters^ It was also manifested in the shape of 
the luminescence line observed in transmission geometry 
with thin samples^ 

In this paper, the minority carrier distribution in the 
samples is obtained by a Monte Carlo modeling of the 
photon recycling process. We assume that the intrinsic 
(unfiltcrcd) emission lincshape is faithfully given by the 
van Roosbroek-Shockley (VRS) relation^ We then pro- 
ceed to calculate the reabsorption probability in terms of 
the measured absorption coefficient and the VRS emis- 
sion spectrum. The only remaining adjustable parameter 
is the quantum efficiency 77 of emission and our experi- 
ments thus provide an accurate measurement of r\. The 
quantum radiative efficiency can also be estimated in- 



dependently from the time-resolved luminescence decay 
experiments^ 7 - and excellent agreement is obtained be- 
tween these two determinations of r\. The index 7 char- 
acterizing the Levy flights distribution is found to be de- 
pendent on the Urbach tailing parameter of the absorp- 
tion spectrum and on the temperature. The observed 
reflection and transmission spectra for differently doped 
InP samples confirm the power-law decay of hole concen- 
tration and thus validate the concept of anomalous hole 
diffusion. 



II. EXPERIMENTAL RESULTS AND THEIR 
INTERPRETATION 

Experimental spectra 

The luminescence spectra were measured in 
moderately-doped (Nd < 10 18 cm -3 ) and heavily- 
doped (up to N D = 8 x 10 18 cm~ 3 ) InP wafers of 
thickness d= 350 [im supplied by Nikko MetalsJ^ For 
some doping levels thinned samples were also studied. 
The spectra were taken at room temperature in both 
reflection and transmission geometries. The excitation 
wavelength A = 640 nm was chosen to ensure short 
penetration of the incident radiation into the wafer, so 
that the resulting distribution of holes is dominated by 
the carrier kinetics. The distribution is slightly modified 
only near the surface in a region whose thickness depends 
mainly on the surface recombination rate. Luminescence 
spectra observed in the reflection geometry for mod- 
erately doped samples with Njj = 2,3 and 5.3 x 10 17 
cm -3 are presented in Fig. [TJ (a). The inset shows the 
same spectra in the logarithmic scale, demonstrating an 
exponential drop-off on both red and blue sides. For 
all moderately doped samples the observed spectra are 
similar, varying only in the luminescence intensity. The 
spectra in both reflection and transmission geometries 
are shown in log scale in Fig. [T] (b) for a sample with 
No = 3 x 10 17 cnr 3 . We have also measured the 
absorption curves in the transparent region up to the 
values a{E) « 200 cm -1 limited by the transparency 
of the samples. To obtain the absorption spectra in a 
broader energy region, we have extrapolated our experi- 
mental data using the available data obtained with much 
thinner samples i 17 ' 20 ' 21 The entire absorption spectrum 
is shown in Fig. [Tj (b) by a dash-dotted line; also shown 
is the intrinsic emission spectrum Svrs, calculated with 
the van Roosbrocck-Shockley relation ) 18 ' 22 

S V rs{E) oc E 3 ai {E)cxp(-E/kT) . (1) 

Here oti(E) is the interband contribution to the absorp- 
tion coefficient a(E) calculated as on{E) = a(E) — a>f c , 
where aj c is the free-carrier absorption coefficient (asso- 
ciated with the interband transitions to the upper con- 
duction valley 2 ^). The ctf c is linear in doping and does 
not tangibly vary with energy in the region of interest. 
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FIG. 1. (Color online) (a) Experimentally observed luminescence spectra in the reflection geometry for differently doped 
samples at T = 296 K. Inset shows the spectrum in the log scale. The right panel (b) shows the luminescence spectra both in 
reflection and transmission geometries for Nn — 3 x 10 17 cm" 3 . The dashed line (VRS) shows the intrinsic emission spectrum 
calculated with the van Roosbroek-Shockley relation, Eq. ([T]). The dash-dotted line shows the absorption coefficient for this 
sample. 



Therefore, its subtraction can be done with sufficient 
accuracy.— 

Figure [T] (b) clearly indicates the red shift of the ob- 
served luminescence spectra relative to the intrinsic emis- 
sion spectrum and the absorption spectrum of the sam- 
ple. Both in the transmission and reflection geometries 
the maxima of observed lines arc on an exponentially 
steep slope of the absorption spectrum. The reflection- 
geometry spectra at low doping (up to 5 x 10 17 cm" 3 ) are 
shifted to the red by about 30 meV as compared to the 
VRS line. Spectra in the transmission geometry show a 
much larger red shift (of « 70 meV) and a steeper decay 
at the high-energy side. The distortion of the intrinsic 
VRS spectrum is due to reabsorption of the outgoing ra- 
diation. 

The effects of reabsorption can also be traced in the 
shifts of the peak energies E max of the photolumincs- 
cence spectra relative to those of the intrinsic spectra. 
These shifts in both reflection and transmission geome- 
tries are illustrated in Fig. [2] for moderately-doped and 
heavily-doped samples. Shown are the i? max values ob- 
tained in this work for differently dopes samples with 
d = 350 fim as well as those observed by Sieg et al£ in 
reflection geometry for epitaxial layers of thickness 2 /xm. 
The epitaxial results exhibit much smaller reabsorption 
red shifts. For comparison, Fig. [2] also shows the E max 
positions for intrinsic spectra calculated with Eq. (fl]). 
The steady increase of E mlix with the majority carrier 
concentration is explained by the conduction band fill- 
ing and the Moss-Burstein shift. The smaller increase of 
_E max at high doping for the transmission spectra is due to 
the smaller spread of the minority carriers, as discussed 
below. 



Interpretation: filtering functions 

The observed spectrum in reflection geometry can be 
described by the following expression 

Irefl(E) = SvRs(E)F re fi(E) , (2) 

which defines the filtering (reabsorption) function 
F re fi (E) . The spectrum observed in transmission geom- 
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FIG. 2. (Color online) Positions of the maxima of the photolu- 
minescence spectra in reflection and transmission geometries 
for differently doped samples. Open dots show the results in 
reflection geometry from epitaxial samples by Sieg et al£ The 
maxima calculated for van Roosbroek-Shockley (VRS) spec- 
tra are shown by full dots. Dashed lines are guides for the 
eye. 
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FIG. 3. (Color online) (a) Ratio of the spectra observed in transmission and reflection geometries, 1Z(E) — Itr(E)/ I re fi(E) 
for moderately doped samples (dots), d = 350 fim. Lines show the ratio TZ(E) calculated for a sample with Nd = 3 x 10 17 
cm -3 , assuming exponential distributions of the minority carriers corresponding to diffusion lengths Ld = 50 /im and 85 fim. 
(b) Ratio 1Z(E) for heavily-doped samples (dots). Solid lines show the results of calculations based on Eqs. (|5I6[) using the 
diffusion approximation for the hole distribution with diffusion length as a fitting parameter, see Table [I] 



ctry can be similarly expressed in the form 

I tr (E) = S VRS (E)F tr (E) , (3) 

where F tr (E) is the filtering (rcabsorption) factor in 
transmission geometry. 

The ratio of experimental spectra, 7Z(E) = 
Itr(E)/ I re fi(E), measured in transmission and reflection 
geometries is shown by the dots in Fig. [3] (a) for three 
moderately doped samples with different doping levels. 
According to Eqs. (|2I3|) . this ratio reduces to the ratio of 
the corresponding filtering functions, 1Z(E) = F tr /F re fi, 
with the factor Svrs canceling out. The ratio 1Z{E) is 
convenient for our analysis becasue it does not depend 
on the details of the iintrinsic VRS spectrum. For all 
samples a steep decline of the 71(E) with E is observed 
in the region from E « 1.27 eV to E w 1.31 eV, with 
a minor shift of the curve to lower energies for higher 
doped samples. 

The absorption and luminescence spectra for heavily 
doped samples were reported earlier Ji. These spectra 
strongly depend on the doping, predominantly due to the 
conduction band filling and the attendant Moss-Burstein 
effect. For heavy doping, the luminescence lineshapes 
are poorly described by the VRS relation. Presumably, 
while in low doped samples the shape of the spectra is 
dominated by the Urbach exponential decay in the red 
wing and the temperature-dependent exponential decay 
in the blue wing, the main input for the heavily-doped 
samples comes from the absorption region a = 10 3 — 
10 4 cm -1 where the behavior a(E) is non-exponential. 
Measurements of the energy dependence in this region 
have a large uncertainty. Experimental results for the 
ratio Ti-{E) for these samples are shown by dots in Fig. 
[3] (b). In this case the ratio variation depends on the 



doping. The Moss-Burstein shift manifests itself in the 
higher-doped samples as a shift in the curves TZ{E) to 
the high-energy side. 

The effect of the Moss-Burstein shift is clearly seen in 
experimental results for the ratio of the total (integrated) 
intensities of the emission lines IZtot = Itot,tr/Itot,refi for 
the whole set of differently doped samples presented in 
Fig. @]as a function of the doping level, TZtot(Nd)- Two 
notable features of this dependence are: (i) high values 
of the transmitted radiation intensity; and (ii) minimum 
of Ktot at N d 2 x 10 18 cm" 3 . The latter effect reflects 
subtle details of the rcabsorption probability (cf. Fig. [5]) 
and will be discussed below. 
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FIG. 4. (Color online) Ratio of the total luminescence inten- 
sities IZtot for differently doped samples, d = 350 ^m. Solid 
line is a guide to the eye. 
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To check the influence on the spectra of the in-plane 
distribution of the excitation, we varied both the size of 
the excitation spot and the size of the registration area of 
the observed spectra. Only minor changes in the spectra 
were observed. This indicates that the only dependence 
of importance is the distance from the surface of excita- 
tion. Hence one-dimensional kinetics model is adequate 
for the spectra interpretation. 

It is important that due to the high refractive index of 
InP, the total internal reflection angle for radiation com- 
ing out into air is 6t — 17.5° [so that cos(9 t ) = 0.954]. As 
a result, the spread of the travel distances for outgoing 
radiation is negligible. We assume therefore that the ob- 
served radiation propagates along the normal to the sam- 
ple surface (we have checked that allowing for the angular 
spread does not lead to sizable corrections) . The filtering 
factors F re fi (E) and F tr (E) can be expressed through the 
one-pass filtering functions F 1 (E) and F 2 (E), defined by 



Fi(E) = f p(z)exp[-a(E)z]dz, 
Jq 



F 2 (E) = / p(z)exp[a(E)(d- z)}dz 



(4) 



where p(z) is the non-equilibrium hole concentration at 
a distance z from the sample surface. 

For a sample of finite thickness, taking into account 
multiple reflections of the luminescence radiation, we 
have 

F refl (E) (1 - R) 1 _ R2exp{ _ 2ad) ■ 



(5) 



and 



F tr (E) = (1 - R) 



F 2 + RFi exp(-ad) 
1 - R 2 exp(-2ad) 



(6) 



Here R is the InP reflection coefficient. For our samples it 
was measured to be R = 0.3 close to the absorption edge; 
in the further numerical calculations we took into account 
the experimental dependence R(E) and its variation with 
doping. 

As seen from Eqs. (|2I3[) and (|5l6p . the ratio 



11(E) 



F tr (E) _F 2 + RFx exp(-ad) 



F reJ i{E) Fi + RF 2 exp(-ad) 



(7) 



has important advantages for the analysis of the minor- 
ity carrier profile. As noted already, this ratio does not 
depend on details of the intrinsic emission spectrum. Fur- 
thermore, it is not sensitive to multiple reflections, since 
the denominators of Eqs. (|5I6[) cancel out. At the same 
time, it is quite sensitive to p(z) through F\(E) and 
F 2 (E). Therefore, the ratio 71(E) can be used to quantify 
the spatial hole distribution. 

For a rough estimation of the hole spread, we have 
calculated the ratio 71(E) assuming that the hole con- 
centration decays exponentially away from the surface, 



TABLE I. Effective diffusion length for differently doped sam- 
ples obtained from the data for thick d — 350 )im samples. 
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40 
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p(z) = p(0) cxp(— z/Ld), where Ld is an "effective dif- 
fusion length" that we use as a fitting parameter. The 
results for moderately doped samples are shown in Fig. 
[3] (a) by dashed lines for two values of Ld- We see that 
Ld = 85 ^m describes the experimental results very well, 
while the smaller value Ld = 50 /im gives a steeper de- 
cay than that observed experimentally. Both values are 
far above the typical hole diffusion lengths (which do not 
exceed 10 /*m). 

For the heavily doped samples, Fig. [3] (b), the ex- 
ponential approximation for p(z) does not fit the exper- 
imental data quite well. The best fit procedure gives 
the values of Ld listed in Table HI The energy varia- 
tion of the ratio is dominated by the strong decay of 
F tr (E) oc exp[— a(E)(d — Ld)], while the decrease of 
F re fi(E) is shifted to higher energies due to the shorter 
escape distance. Therefore, any distribution with an av- 
erage spread length Ld that is considerably shorter than 
the layer thickness d will produce curves similar to those 
depicted in Figs. [3] As seen from Tablc[H the hole spread 
in the layer is reduced in heavily doped samples but still 
is much broader than could be anticipated for a conven- 
tional hole diffusion, indicating the importance of recy- 
cling effects. 

Note that the exponential decay approximation with 
an enhanced effective diffusion length fails to describe the 
variation of 71(E) with the sample thickness, especially 
for thin (« 50 /jm) samples. 

As a matter of principle, accurate measurements of 
the reflection and transmission spectra together with the 
absorption spectra enable us to find F re fi and F tr and 
from these functions recover p(z). Indeed, using the 
experimentally known function a(E) one can calculate 
Svrs(E). As discussed above, this step works well for 
moderately doped samples. Next, we find F re fi(E(a)) 
and F tr (E(a)) via Eqs. (|2I3[) . From these functions one 
can find the one-pass filtering functions (g]), regarding 
them as functions of a, viz. 



Fi(a) 
F 2 (a) 



F re fi - RF tr exp(-ad) 
1 - R 

F tr - RF ref i exp(-arf) 
1 - R 



(8) 



Equations (JSJ) with their right-hand sides expressed in 
terms of the experimental data represent experimental 
determination of the one-pass filtering functions. Then 
Eqs. (|4|) can be viewed as integral equations for p(z). 
In fact, they represent a Laplace transformation (with 
the argument a) of the hole distribution function p r (z), 
defined (for F\) as p r (z) = p(z) for < z < d and 
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FIG. 5. (Color online) Reabsorption probability V(z) for moderately doped samples (a) and for heavily doped samples (b) 
calculated with Eq. ([9]) (dots). Dashed lines correspond to the power-law approximation, Eq. (|12[) . 



p r {z) = otherwise. Similarly, for F2 one can consider 
Pt{z) = p(d — z) inside the sample. In principle, one can 
obtain p(z) by inverting the Laplace transformations nu- 
merically, using one of the developed approaches^ The 
pair of Eqs. ([8]) can be used for a subsequent reiteration 
scheme in the numerical solution. However, the numer- 
ical inversion of the Laplace transformation is a classi- 
cal ill-conditioned mathematical problem^ Attempts to 
find p(z) in this way suffer from a low solution accu- 
racy. The dominating errors in the filtering functions 
appear in the red-wing region of the experimental spec- 
tra, where the luminescence intensity is low and the noise 
signal is high. This high-noise signal in the low-energy 
region (where Fi « F2 « 1) is clearly seen in Fig. [3] 

Therefore, instead of attempting to invert the Laplace 
transformations, we have modeled the distribution of 
holes in the layer, calculated the filtering functions, and 
compared the results with the functions found from the 
experiments. 



TABLE II. Parameters of the reabsorption probability distri- 
bution ||12p and the recycling factor $ for differently doped 
samples; the values are from Ref. Il7l 



No, 


10 18 cm" 3 


0.2-0.6 
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3.7 


6.3 


8 


7 




0.64 


0.79 


0.7 


0.64 


0.69 




, ^m 


0.0962 


0.625 


0.7 


1 


1.377 


$ 




90-97 


32 


IS 


10 
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#1 






48 


20 


12 
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Here Af(E) is the normalized spectral density for the 
number of photons in intrinsic luminescence, J\f(E) = 
CSvrs(E)/E, cf. Eq. ([T]), and C is a normalization 
constant. The exponential integral function Ei(l, \z\) is 
defined by 



, . f°° dt 

Ei (M) = y y exp 



(-**) 



(10) 



III. MINORITY CARRIER DISTRIBUTION IN 
A FINITE LAYER 



Studies of the temporal relaxation of non-equilibrium 
minority-carrier excitations in doped rt-InP demonstrate 
a very high quantum radiation efficiency and high val- 
ues of the photon recycling factor— (which in the bulk 
crystal equals the ratio of radiative to non-radiative re- 
combination rates). Therefore the hole distribution in 
case of the stationary excitation is strongly influenced by 
the recycling process. In an infinite medium with the 
initial excitation homogeneously distributed in the plane 
z = z' the probability of reabsorption in the plane z = z" 
is given by&l^ 



V{z"-z') = l - f Af(E)ai(E)Ei[l,a(E)\i 



z'\]dE. (9) 



Due to the residual absorption the full probability of in- 
terband reabsorption P to t < 1- However, in the range 
of thicknesses and distances of interest the difference is 
negligible for all our samples, so that 



V(z)dz 



'm)^E 

a(E) 



1 



(11) 



Results of numerical evaluations of "P(|z|) with Eq. (|9]) 
are shown in Fig. [S] (a). For the entire range of z, the 
dependence "P(|z|) is very close to 



V{\z\) 
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ryz ■ 
1 mm 



2(z„ 



+1 



(12) 



For moderately doped samples, z m i n ~ 0.1 yum and the 
index 7 sa 0.65 while slightly decreasing with the doping 
level. 
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FIG. 6. (Color online) (a) Stationary hole distribution p(z) (in log-linear scale) calculated by Monte Carlo with the reabsorption 
probability (|12[) for Nd = 6.3 x 10 18 cm~ 3 and different boundary conditions at the back surface (r: reflective, a: absorptive, 
and nr: non-reflective). Dashed lines show the fitting with the power law p(z) = c/(z + z m m) 1+7 . (b) The hole distribution 
(in log-log scale) for differently doped samples; the sample parameters are listed in Table HT1 The dashed line corresponds to 
an exponential distribution with Ld = 45 /im (cf. Table [J). 



The emission-reabsorption process described by Eq. 
(J9j) can be interpreted as a one-dimensional random walk 
of holes. Its stochastic nature is associated with the emis- 
sion spectrum M[E). The flight probability, accurately 
described by Eq. (|12[) , is typical for the "anomalous diffu- 
sion" of the Levy-flights type. Its hallmark is the asymp- 
totic spatial decay Viz) oc z" 1 " 7 with a "heavy tail", 
i.e. the power-law asymptotics with an index 7 < 2. 
Moments of this distribution diverge, and so does the 
conventionally defined diffusion coefficient. 

Analytical calculation of the index is possible in a 
model where the interband absorption spectrum is ap- 
proximated by a simple function 

^ = 7T WW FWA1 ' (13) 

1 + cxp[(£ g - E)/A\ 

which decays exponentially below the absorption edge 
and saturates above it. This model accounts for the Ur- 
bach tailing (described by the tailing parameter A), but 
it does not describe the growth of ai(E) at E > E g . 
Similarly, the emission spectrum can be approximated 
by the VRS relation Q [with a t {E) given by Eq. (|I3]t] 
in which the pre-exponcntial factor for A < kT <C E g 
can be replaced by its value at E g . In this model one 
obtains 



Equation (fT4|) predicts a decrease of 7 at lower temper- 
atures. It also explains the decrease of 7 with increasing 
Nd . The later effect is due to the smearing of the absorp- 
tion edge (described by an increase of the Urbach tailing 
parameter A). Estimation of 7 with Eq. (fT4|) for mod- 
erately doped samples (A = 9.4 meV)^ gives 7 ps 0.64, 
close to the results obtained with more accurate numer- 
ical calculations. 



For heavily doped samples, the distributions P(|z|) cal- 
culated numerically with Eq. @ arc shown in Fig. [3] (b) 
with the distribution parameters listed in Table [TTJ The 
calculated variation of 7 is more complicated than that 
predicted by Eq. (|T4|) due to the competition of several 
factors. Thus, compared to the low-doped samples, the 
index increases to 7 = 0.78 for Nd = 2 x 10 18 cm" 3 and 
then decreases for higher doping. For Nd = 6.3 x and 
8xl0 18 cm -3 , the combination of a strong smearing of 
the absorption edge and the Moss-Burstein shift reduces 
the index 7 to 0.64. Besides, the normalized reabsorption 
probability is redistributed to higher distances, leading to 
a broadened hole distributions in the layer both at low 
and high ends of the concentration range. This corre- 
lates with the luminescence ratio dependence shown in 
Fig. [4J namely with the observed minimum in lumines- 
cence intensity in transmission geometry for samples with 
Nd = 2 x 10 18 cm" 3 . Wc remark that the accuracy of the 
■p(|z|) calculation is not high for heavily-doped samples 
due to the less reliable approximation for the emission 
spectrum. 

To calculate the stationary hole distribution, consider 
the temporary evolution of the distribution p(z, t) after a 
short excitation pulse at t — 0. Given the p(z, t), we can 
evaluate the stationary distribution p s t(z) for constant 
excitation rate using the Duhamel principle^ viz. 

/■OO 

Pst (z)= p(z,t)eM-t/r)dt/T , (15) 
Jo 

where r is the average lifetime of holes. 

Equation (fT5|) suggests an approach to the Monte Carlo 
calculation of p s t(z). First, we simulate the distribution 
p(z, N) after a fixed number TV of recycling events of holes 
initially excited at the surface. This distribution can be 
viewed as a distribution p(z,t) of the hole spread for a 
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FIG. 7. (Color online) Spectral ratio 1Z(E) for samples of different thickness d. (a) Results for moderately doped (No = 2 x 10 18 
cm -3 ) samples with d = 50 and 350 /an (solid curves) are compared against model calculations with $ = 32 (dashed lines), (b) 
Results (solid curves) for heavily doped (Nd = 6.3 x 10 18 cm -3 ) samples with d = 250 and 350 /an are fit by model calculations 
(dashed lines) with $ = 10. Back side is assumed to be non-reflecting. 



given time, t = Nr r , where r r is the radiative emission 
time. 

Next, we average the Monte-Carlo simulated distribu- 
tions p(z,t) over the random walk times t distributed 
as exp(— t/r)/r. This average is equivalent to averag- 
ing over the number of recycling events with the average 
recycling factor $, 

oo 

Pst {z) = £ p{z, N) exp(-iV/$) . (16) 

iV=l 

For $ ^> 1 the main contribution to the sum at all z > 
z min comes from the terms with N 3> 1, so that Eqs. 
([TBI and (fl"5|) give very close results. 

The simulated stationary hole distributions are sen- 
sitive to the boundary conditions. To discuss this ef- 
fect, we consider the results obtained for a sample with 
N D = 6.3 x 10 18 cm" 3 , see Fig. © (a). The normal- 
ized distributions are found to be practically insensitive 
to the boundary conditions on the front surface. Re- 
sults presented in Fig. |H] are calculated with the re- 
flective front (rf) boundary condition. In this case the 
stationary distributions satisfy p r f(z) ~ 2p OQ (z), where 
Poo(z) is the distribution for a localized source in an in- 
finite medium. This corresponds to the contribution of 
an "image" source of photons provided by the reflecting 
boundary. 

In contrast, the distribution is substantially affected by 
the conditions at the back surface, see the curves 1-4 in 
Fig. [5] (a). Depicted are the hole distributions for: reflec- 
tive (r), absorbing (a), and "non-reflective" (nr) bound- 
ary conditions (the latter corresponds to a semi-infinite 
sample without a back surface). All distributions show 
a characteristic non-exponential drop-off with distance 
from the front surface, typical for Levy flights. 



In case of reflective back surface, the decline of concen- 
tration appeared to be much slower, p r \, oc z 1+1 where 
7' < 0, as compared to that for the semi- infinite medium 
(or non-reflective back surface, curve 4). The slower de- 
clining distribution is in contradiction with the fast de- 
cline of Fi (E) in the region of a above 200 cm -1 and also 
with the drop-off of 1Z(E) observed in the 50 /jm sample. 

For cither the absorbing back surface or the non- 
rcflcctivc boundary condition, the distribution has power 
asymptotics p(z) tx z" 1 ^ with 7 ps 0.12 < 7. It still ex- 
tends to much wider region than an exponential distribu- 
tion and decays slower than the reabsorption probability 
V(z). The assumption of non- reflective back boundary 
gives a much closer fit to out experimental data for the 
filtering functions in all samples. We have no physical 
reason to assume a totally light-absorbing back surface. 
But it seems reasonable that the surface scatters part of 
the luminescent radiation diffusively, so that it escapes 
from the observation area and does not contribute to fur- 
ther recycling process. The luminescence escape from 
the observation area manifests itself also in the filtering 
function dependence in the transparency region at the 
red wing of the emission line. Therefore, in our compari- 
son with experimental data for the filtering functions and 
their ratios, we shall use the distributions obtained with 
"non-reflective back boundary" . 

The simulated distributions p{z) for the set of differ- 
ently doped samples, assuming non-reflective boundary 
conditions at the back surface, are shown in Fig. [S] (b) in 
log-log scale. For comparison, we also show an exponen- 
tially decaying distribution corresponding to Ld = 45 fjja 
(cf. Table |J). We note that the non-exponential decay of 
the concentration is more sensitive to the recycling factor 
than to the variation of 7. 

Using the hole distributions of Fig. [5] and Eqs. d5j CEO) ^ 
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we calculated TZ(E) for all samples, using $ as the only 
adjustable parameter. The best-fit values of $ are also 
listed in Table HH These values are close to but some- 
what smaller than the values obtained independently 
from the time-resolved luminescence studies^ The most 
sensitive comparison of the theory and experiment is pro- 
vided by the samples with the same doping Nd = 2 x 10 18 
cm -3 but very different thickness, see Fig. [7] (a). Cal- 
culations with the modeled anomalous distributions give 
an excellent agreement with the experimental data. Even 
though for the relatively thick samples [d=250 and 350 
/jm, see Fig. [7J (b)] a satisfactory agreement can be also 
obtained assuming an exponential hole distribution, this 
assumption fails completely when one considers a pair of 
samples of vastly different thickness (e?=50 and 350 fim). 
The exponential approximation fails to describe thin and 
thick samples simultaneously. 

The low-doped samples permit a more detailed com- 
parison of the experimental and calculated spectra. If we 
assume that the initial spectra are well described by the 
van Roosbroeck-Shockley relation (TTJ), then we can use 
Eqs. (|5I6[) to find the filtering functions for the trans- 
mission and reflection spectra separately (as well as their 
ratio). Comparison of the filtering functions for the sam- 
ple with n — 3 x 10 17 cm -3 is shown in Fig. [5] The 
results for the ratio do not differ noticeably from the ex- 
perimental data shown in Fig. [3j In the rapid decay 
region of the filtering functions, the agreement with ex- 
periment is quite good. However, there is a deviation 
in the region of small values of a(E). One should note 
that in this region the variation of functions F tr (E) and 
F re fi (E) is sensitive to multiple reflections of the outgo- 
Absorption coefficient, cm 1 
1 10 100 1000 10000 
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FIG. 8. (Color online) Filtering functions F tr (E) and 
F re fi(E) (blue and red dots, respectively) for Nd — 3 x 10 17 
cm" 3 sample and their comparison to model calculations. 
Hole distribution is obtained by Monte Carlo modeling with 
the reabsorption probability given by Eq. (|12[) and the quan- 
tum efficiency r\ = 98%, ignoring multiple reflections; dashed 
lines correspond to allowing multiple reflections. 



ing radiation from both the front and the back surfaces. 
Again, some of the angle-distributed radiation goes out 
of the observation area and is therefore lost. In the ex- 
treme case of a narrow observation area, only one-time 
reflection from the back surface contributes to the spec- 
tra. The filtering functions calculated barring multiple 
reflections are shown in Fig. [5] by the solid lines, while 
calculations allowing all multiple reflections are shown 
by the dashed lines. Since the experimental data in the 
red wing of the dependences F tr (E) and F re fi(E) fall 
between these lines, we conclude that the discrepancy is 
likely due to a deviation of the experimental conditions 
from the assumed one-dimensional geometry. This, of 
course, does not invalidate the results in the wider en- 
ergy range, where the filtering functions exhibit rapid 
variation. 



IV. DISCUSSION 

As was shown above, the experimental results for the 
filtering functions and their ratio are well interpreted in 
a model where the only mechanism of the hole spread 
in the sample is by photon recycling. Effects associated 
with the conventional mobility of holes reside in the pa- 
rameters of the model, such as the recycling probabil- 
ity distribution. Thus, the minimum recycling distance 
incorporates the hole diffusion length and the recycling 
factor $ reflects the losses through all channels including 
surface recombination. 

The main result is that the hole concentration does not 
decay exponentially from the wafer surface, as it would 
for a normal diffusion process. The spread resulting from 
the photon-assisted random walk of holes is in agreement 
with the Levy-flights theory, which predicts heavy power- 
law tails in the distribution. 

It is interesting to compare the results for the quan- 
tum radiative efficiency 77 and the corresponding recy- 
cling factor obtained in this work, with those obtained 
earlier— from the kinetics of time-resolved luminescence. 
The concentration dependence of these very important 
parameters is illustrated in Fig. [9] for both the earlier 
and the present studies. Both sets of results are also 
listed in Table |TT] in the instance of the recycling factor, 
as and respectively. The earlier data were obtained 
with the excitation photon energy in the relatively low- 
absorption part of the spectrum (red wing of the funda- 
mental absorption edge) and therefore they correspond 
to a much larger penetration depth of the excitation and 
a different resulting hole distribution. Besides, in the 
earlier experiments^ the escape of photons from the ob- 
servation spot on the sample surface was found to be of 
importance, providing and additional mechanism of the 
photon loss. 

The data for the quantum radiative efficiency 77 from 
Ref. [I?] are shown by the circles in Fig. [T3 The solid line 
drawn through these circles is an interpolation based on 
the standard decomposition of hole recombination rates 
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into the radiative (B) and non-radiative terms, impurity 
(A) and Auger (C), viz. rj(N D ) = BN D /(A + BN D + 
CNp). The solid line for $ shows the same data re- 
calculated as $(Nd) = 77/(1 — 77), which gives an upper 
estimate for the recycling factor M- The squares show the 
results of the present work. Both the high values of the 
recycling factor and its concentration dependence are in 
close agreement with the earlier results, obtained by a 
very different technique. 

The most direct experimental manifestation of the 
non-exponential decay of the hole concentration is found 
in the ratio TZ(E) of the spectra observed in transmis- 
sion and reflection geometries for the samples of the same 
doping level but very different thicknesses, see Fig. [7] In- 
deed, the ratio Tt{E) will rapidly decline in energy only if 
the distribution itself decays within the thickness of the 
sample, as is supported by the data for 50 ^m-thick sam- 
ples, which suggest that the spread size is shorter than 
the sample thickness. This, however, is in contradiction 
with the fairly high intensity of transmitted radiation 
for much thicker (350 /im) layers of the same doping. 
The power-law decay of the hole concentration is steep 
enough at short distances (steeper than an exponent) to 
fit the data for the thin sample, and at the same time 
slow enough at large distances (again, compared to an 
exponent) to account for the data for thick samples. The 
set of samples with low doping but different thicknesses 
should be most useful for further detailed investigations 
of the anomalous transport. 

Our experimental data can be well described by the 
distributions obtained with different boundary conditions 
at the front (excitation) surface, since the concentration 
decay is relatively insensitive to these conditions. This 
is not the case for the back-surface boundary conditions. 
Our results arc not in agreement with the assumption 
of almost total reflection from the back surface that one 
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FIG. 9. (Color online) Quantum radiative efficiency r\ (circles 
and solid line) from Ref. [l?] and the recycling factor $, both 
calculated from the earlier data (solid line) and obtained from 
the present experimental data (squares), cf. Table UTI 



could expect based on the refractive index ratio. We at- 
tribute this discrepancy to the radiation escape from the 
observation region, in agreement with our earlier con- 
clusions in time-resolved photoluminescence studies^ 
Good agreement is obtained for either absorbing or "non- 
reflecting" back surface. 

Our transport problem is well approximated by a sta- 
tistical model that has been thoroughly studied, namely 
that of a one-dimensional random walk on an infinite 
lme&2& Therefore, apart from some modifications at the 
back surface (see Fig. , a suitable distribution of holes 
after N S> 1 flights originating from z = should ap- 
proach the so-called standard stable distribution^ with 
zero asymmetry and a given index 7, 



p(z,N) 



cos{kz)e- N{z " kr ' dk 



(17) 



where z c is the depth scaling factor. By the order of mag- 
nitude one has z c = z m i n (comparison of the distribution 
(fTT|) with our Monte-Carlo results for moderately doped 
samples gives z c = 0.23 /mi). To find the stationary dis- 
tribution for a given finite average number $ of random 
flights, one should use Eq. (fT5"j) and integrate Eq. (fTT|) . 
resulting in 



p(z,F) 



cos(kz) 

1 + Wzjcyi 



dk . 



(18) 



Equation (fT5)) yields a distribution that is close to our 
Monte Carlo results. Furthermore, it can be analyzed 
quantitatively. Analysis of Eq. (fT8]l readily shows the 
existence of two regions in the excitation front spread, 
that of asymptotic decay and that of short flights. In the 
asymptotic region, the hole concentration is dominated 
by the one-pass long-distance flights and hence decays as 
the reabsorption probability itself, p(z) oc z _7_1 . The 
asymptotic region corresponds to z > z^ 1 ! 1 = Zf. For 
$ ;§> 1 and 7 < 1 the front spread distance Zf ^> z c . 
In the region z <C Zf, large values of k contribute to 
the integral in (|18[) . Then one can neglect unity in the 
denominator of (|18|) . which gives p(z) oc z 7_1 with much 
slower decay with distance than in the asymptotic region. 

For moderately doped samples, one has z/ ~ 300 /im 
and the onset of the asymptotic dependence is at dis- 
tances larger than the sample thickness. Experimentally, 
neither d -C Zf nor g? 3> Zf conditions apply. This ham- 
pers the determination of the index of anomalous dif- 
fusion through the luminescence spectra. The experi- 
ments for the radiation escape from the remote edge of 
the wafer— are better suited for this goal. 

Finally, we note that the one-dimensional Levy flights 
transport is described by an integro-differential equation, 
corresponding to the diffusion equation with the recy- 
cling termi^ This equation can be solved numerically, us- 
ing available COMSOL software. Evaluation of the hole 
distribution in this way gave excellent agreement with 
the Monte Carlo results, except within an area of the 
order of the hole diffusion length near the sample sur- 
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face. In the far-distance asymptotic region, the integro- 
differential equation admits of an analytic solution^ 
which is again the distribution (|18|) . 

V. CONCLUSION 

The luminescence spectra of moderately thick n-type 
InP wafers have been studied in reflection and transmis- 
sion geometries. Analysis of the reflection and transmis- 
sion spectra and especially of their ratios provides an 
unambiguous evidence for the anomalous non-diffusive 
transport of minority carriers via photon recycling. As 
a result, the minority-carrier distribution acquires heavy 
tails, accompanied by an anomalously high intensity of 
luminescence in the transmission geometry Our inter- 
pretation of the experiment provides an independent es- 
timate of the photon recycling factor $, which is in good 
agreement with earlier results obtained by time-resolved 



luminescence studies. 

The analysis also confirms the very high quantum radi- 
ation efficiency in moderately doped samples. Such sam- 
ples provide, therefore, a remarkable experimental sys- 
tem for studying the Levy- flight transport. Its unique 
advantage is the possibility to vary the Levy distribution 
index by changing the temperature or the doping. 
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